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Abstract 

The dynamics of preheating after inflation has not been clearly understood yet. In 
particular, the issue of the generation of metric perturbations during preheating on 
super-horizon scale is still unsettled. Large scale perturbations may leave an imprint 
on the cosmic microwave background, or may become seeds for generation of primor- 
dial black holes. Hence, in order to make a connection between the particle physics 
models and cosmological observations, understanding the evolution of super-Hubble 
scale perturbations during preheating is important. Here, we propose an alternative 
treatment to handle this issue based on the so-called separate universe approach, 
which suggests less efficient amplification of super-Hubble modes during preheating 
than was expected before. We also point out an important issue which may have 
been overlooked in previous treatments. 

1 Introduction 

Rather recently, people realized that the process of reheating after inflation need not simply be described 
by the usual single particle decay of the inflaton field. If there are relevant scalar field degrees of freedom 
coupled to the inflaton field at the end of inflation, the coherent oscillation of the inflaton can excite 
those scalar fields by parametric resonance 1 . This process was extensively investigated, in mast cases, 
neglecting metric perturbations, which typically will be a good approximation as long as we are interested 
in processes occurring within the Hubble horizon scale. Often the energy transfer due to parametric 
resonance is much more efficient than that due to the usual particle decay, and sometimes it is violent. 
It was then suspected that there might be significant effects on the large scale perturbations beyond the 
Hubble scale[21 due to resonance. Some particle physics models may predict distinguishable signals in 
the cosmic microwave background, or formation of inconsistently large number of primordial black holes 

m 

Here in this short contribution, we would like to propose another approach to discuss this issue 
of the generation of large scale perturbations at preheating. We make use of the separate universe 
approach We stress that unphysical acausal energy transfer, which may not carefully prohibited in 
previous treatments, does not take place in this new approach. 



2 Separate universe approach 

Basic Idea : We consider the system composed of many scalar fields. The discussions in this paper is 
based on the following statement obtained in Ref. |S10J|2]. 

When we solve the background spatially homogeneous field equation, the e-folding number N = 
log a until the expansion rate H reaches a given value H r depends on the initial value of the fields 
<fr (and their time derivatives). Let's denote this value of N by Nu r - Suppose that the initial value 
of €> has spatial variation over super-Hubble scales. Then, the value of Njj r becomes a function of 
the spatial position x. In this case, the newly generated curvature perturbations on the constant 
Hubble slice TZ is simply given by AA^, i.e., the fluctuations in Nn r . 
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Mechanism of Generating Curvature Perturbations from Isocurvature Perturbations : 

Let's consider matter fields consisting of two components. We call them <j> and \. They are not 
necessarily scalar fields here. We assume that the rates of change of the energy density of <f> and x due 
to the expansion of the universe are different. Namely, 

dlogP0 Qn . v d\ogp x 
dN =-3(l + W0), m = -3(1 + w x ), (1) 

For example, the energy density of <f> changes like dust matter, while x behaves as a radiation fluid. In 
this case, we have — and w x — 1/3. Let us denote the fraction of the energy density of x by 
/ := p x / p. Then, we have 

^ = -3(l + ^) + 3(^- % )/. (2) 

From this expression, we notice that there are two factors necessary for generation of super-Hubble 
adiabatic perturbations. One is non- vanishing w^, — w x , and the other is inhomogeneity of /. 

3 Conformal two field model 

Model : We consider the model whose potential is given by 

V = ^+ 1 - 9 W. (3) 

When g 2 = 2A, the longest wavelength mode of x is the fastest growing mode, and the initial spectrum 
of x is scale- invariant. This appears to lead to exponential growth of pre-existing entropy perturbations 
even when backreaction is taken into account j7J . The equations of motion for the spatially homogeneous 
background are given by 

& + 3H<S> l - — = Q, H 2 = — (-^ 2 + v). (4) 

Here we used the vector notation to represent the two component scalar fields, <p and x, simultaneously. 
Namely, $ x = 4> an d $ 2 = x- 

Simplified Analysis with Time Averaging : We define the averaged energy in a unit comoving 

volume by £ = ^a 3 (h^ + V^j \ , where (• ■ •) represents time averaging. Let us assume that «& 2 is 

bounded and is oscillating. Then, we have many zeros of If we take the average between two zeros, 

we have 

Since H is slowly changing, we obtain an approximate relation: (^a 3 — A</> 4 — 2g 2 x 2 4> 2 ^j ^ = 0. There- 
fore we have 



From this fact, £ w e N is deduced, and hence we have H 2 oc p w £/a 3 cx a 4 . Anyway, Njj r is basically 
independent of the initial value of the field. Therefore it seems that additional curvature perturbations 
are not generated when the wavelength of a mode is longer than the horizon scale. In the language of the 
preceding section, uu — w x vanishes approximately in this model. Hence, amplification of fluctuations of 
O(l) is hardly expected. Notice that H 2 is roughly proportional to $ 4 but the mass squared is to $ 2 . 
Hence, as the amplitude of $ decays, the mass squared tends to be larger than H 2 . 

Detailed Analysis without Time Averaging : We solve the above set of equations numerically 
varying the initial value of x- Then, we can evaluate the initial value dependence of Nu r , the e-folding 
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number until the total energy density reaches a certain specified value. However, the original model has 
the following shortcoming. As there are only two degrees of freedom corresponding to the homogeneous 
modes of ip and x-fields, the exchange of energy between them continues for ever. Then, occasionally 
the value of <p becomes exceptionally small. In this phase, the period between two neighboring zeros of 
<!><I> becomes very long. If the specified target energy density is reached during this phase, AiV becomes 
largely dependent on the choice of the final value of the Hubble rate H r . Then, interpretation of the 
results becomes difficult. In more realistic situations, the energy density carried by these long-wavelength 
modes will be redistributed to a large number of other degrees of freedom (shorter wavelength modes of 
these or other fields) sooner or later through rescattering. 

Here we assume two radiation fluids, and r x . decays into the component with decay constant 
Tj. The modified set of equations are given by 

(TV K 2 t 1 9 . 2 \ 

¥ + (3H + Ti)** -^1=0, p ri = ±Hp ri + IM^) 2 , H 2 = - U# + V + £>, J ' ( 7 ) 

where p Ti is the energy density of the i-th component of radiation fluids. 

We consider the case with g 2 — 2A = 0.02 and = T x = 10~ 5 . We set the initial condition at 
4> = 1 with <fi — x = 0. Using the convention k 2 /3 = 1, the initial Hubble rate is H = 0.05, where 
we assume a negligible contribution to the initial energy density from the x-field. We integrate Eqs. J7|) 
until the Hubble rate goes down to 10~ 6 , at which point the energy is almost completely transferred 
to the radiation fluids. In Fig. 1, we give a plot of the e-folding number N between H = 0.05 and 
H = 10~ 6 varying the initial value of x- We see small fluctuations of the e-folding number N depending 
on the initial value of %. Each plot contains only 100 sampling points with equal horizontal spacing. The 
amplitude of fluctuation is not 0(1) but still it is significant compared to that observed in the cosmic 
microwave background. An interesting and important point is that the dependence on the initial value 
of x °f this fluctuation is almost random. Two panels in Fig. 1 are plots with different range of the ini- 
tial value of x- Although the scale of the horizontal axis is very different, these two plots look quite similar. 




X,*W 5 iogioZ- 

Fig. 2: Similar plots to Fig.l but the vertical axis is 
Fig.l: Plots showing the initial value dependence of the fraction of the radiation energy originating from 

the e-folding number N. In each panel 100 discrete the secondary field \- The panel (b) is a close-up view 

points are plotted. The scale of the horizontal axis is of the panel (a) for small initial value of X- We can see 

different, but these two plots look quite similar. This a trend that the fraction increases as the initial value 

will mean that the distribution of N is rather random. of \ increases. 

The origin of these fluctuations can be understood as follows. If the oscillation of fields is fast compared 
to the change rate of H, the approximation using the time averaging must be good. The error will be 
exponentially suppressed. As mentioned above, however, the period of oscillation occasionally becomes 
very long if the motion of the scalar fields is trapped by chance in the region where the amplitude of 
oscillation in the direction of 0-field is exceptionally small. In such situation, the approximation we 
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made before is no longer valid. Since the motion of the two scalar fields are rather chaotic, the initially 
neighboring trajectories in the configuration space of <j) and x soon deviate. Thus, it is quite sensitive 
to the tiny difference in the initial condition whether this entrapment occurs or not. As a result the 
dependence on the initial condition is almost completely randomized. 

Then, the question is whether this amplification of fluctuations can be a process that generates 
significant fluctuations at the super-Hubble scale during preheating? Since the amplitude of the generated 
curvature perturbation is so sensitive to Xii the variance seen in Fig. f will basically represent the variance 
of random fluctuations incoherent beyond the Hubble scale. If the fluctuations are completely random 
beyond a certain cutoff scale A:" 1 , we have n = 4 spectrum for k < k c in general. (Here we use the 
notation that n — 1 corresponds to the scale invariant spectrum.) 

However, it is not certain whether a smoothing of N(x), e.g., N(Xi) '■= J dxN(x)W(x — Xi) with an 
appropriate window function W is tilted or not. Hence, the averaged value on very large scale may have 
an additional fluctuation other than the n — 4 random component. The amplitude of this additional 
component is difficult to determine because we need to reduce the contamination of n — 4 component. 
To do so, we need to compute the average over sufficiently large volume. This means that we need to 
calculate the value of N for sufficiently large number of samples changing the initial value of x, which 
we have not done yet. Hence, we have not completely rejected the possibility of generating fluctuation 
coherent beyond the Hubble scale. However, the amplitude of this coherent fluctuation should be much 
more suppressed than the incoherent one. 

Generation of Curvature Perturbations -. 

If the exact conformal invariance is broken by some reason, we can expect a curvature perturbation 
may be generated. The easiest way will be changing the decay process of the two fields. In the above 
calculation, we have introduced the decay terms by hand assuming that the both scalar fields decay into 
a radiation fluid with the same decay constant T. We distinguished two components of the radiation 
fluids depending on which scalar field becomes the origin. Suppose, one radiation fluid is composed of 
massive fields, although the mass is negligible during the preheating epoch. Then, sooner or later, the 
mass of this radiation fluid becomes larger than the Hubble rate, and the energy density of this field 
starts to behave like cx e~ 3N , while the other continues to evolve as oc e~ AN . Therefore the isocurvature 
perturbation imprinted as the ratio £ :— Pr x /{Pr 4 , + Pr x ), is transformed into the curvature perturbation 
at this stage. 

Now let us look at the plot of this ratio as a function of the initial value of x- From Fig. (2a), it seems 
that the dependence on Xi is just random fluctuations. However, if we close up the portion for small 
\Xi\ (Fig. (2b)), we can see clear dependence on Xi, which can be explained as follows. For small Xii the 
decay to the radiation fluid becomes effective before the equi-partition between the homogeneous parts 
of the </> and x fields is achieved by the parametric resonance. In this case, the decay to the radiation 
fluid is more efficiently through the <f> field. As a result, the ratio £ becomes small for small \xi\- This 
mechanism to generate large scale fluctuations is very efficient because the fluctuation continues to be 
amplified until the matter radiation equality. Hence, this process will be a possible mechanism to create 
primordial black holes. 

4 Discussion 

In this paper, applying the separate universe approach, we discussed generation of metric perturbations 
at super-Hubble scale during preheating. As an example, we considered the conformal two scalar model 
whose potential is given by Eq. (PJ. Previously, it was argued that large amplitude of fluctuations can 
be generated during the reheating process, but we found that the separate universe approach suggests 
relatively suppressed amplitude of fluctuations. 

We suspect that the origin of this discrepancy is in an overlooked point in previous analyses. Usually 
people decompose fields as (homogeneous background)-!- (perturbations), using the common homogeneous 
background for separated universes. Let's consider two largely-separated homogeneous patches, universe 
A and universe B. Suppose that the preheating is efficient in the universe A but inefficient in the universe 
B. Then, we physically expect that the back reaction to the inflaton oscillation is also efficient in the 
universe A. However, as long as we use the common homogeneous background, the back reaction works 
in the same way both in the universe A and in the universe B. This means that artificial and acausal 
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energy transfer from the universe B to A easily takes place in the treatment using the homogeneous 
background. 

However, the above criticism to the previous works does not apply to the numerical studies. For 
example, in Ref. |SJ, it was argued that super-Hubble modes are enhanced even in a single field model. 
They considered X(j) 4 model. When they simulate the case that only a single mode in the resonance band 
is excited, there was no enhancement in super-Hubble modes. On the other hand, when all the modes 
in the same lowest resonance band are excited, there was enhancement in super-Hubble modes. Their 
explanation is as follows. The gauge invariant perturbations at super-Hubble scale are generated by the 
quadratic product of smaller scale perturbations in the resonance band. Namely, Ck oc J d 3 fcx(k')x(k— k'). 
This integral becomes approximately constant for small k less than the width of the resonance band Afc. 
For larger fc, this integral is small because, when one of the arguments of x is m the resonance band, 
the other is outside of it. So this enhancement mechanism does not apply for small scale modes. In this 
sense, the large scale modes are selectively enhanced. 

However, since £k cx fc -3 / 2 for the scale invariant spectrum, constant value of (k does not mean 
enhancement of larger wavelength mode. Therefore the mechanism of the enhancement of super-Hubble 
modes works only when the cutoff scale aAfc -1 is much larger than the Hubble radius. We can easily see 
that this condition is not satisfied in this model. For the redefined field = a<fi, the field equation becomes 
4>" + X(f) 3 — (a"/a)<f> — 0. Here we used the conformal time coordinate rj = Ja" 1 ^, and a prime denotes 
differentiation with respect to r\. Since the radiation dominant universe is a good approximation for this 
model, we use a = (t/to) 1 / 2 . Then, 77 = 2i a = 1/H = 1/aH, and therefore H = 2to/ri 2 . The amplitude 
of oscillation of <f> is constant. We can set the amplitude of A</> 2 to unity without loss of generality. In 
fact, using the Friedmann equation, we have A</> 2 « (87rG/3a 2 )A0f = TC 2 — l/(2toa) 2 , where we have 
used the fact that the initial value of 0, <P,,. is 0(m p i). Choosing to appropriately, we can set (], = 1. This 
also means that the comoving wave number corresponding to the initial Hubble horizon Ti, is 1. 

The resonance effect becomes important after fj,i] becomes O(l), where fj, is the growth rate of the am- 
plitude due to parametric resonance. Hence, 77 ~ /i -1 when the parametric resonance becomes important. 
At this epoch, the comoving horizon scale is TC = r/ _1 w /j. Then, the condition for the enhancement of 
the super-Hubble modes will be \x ^> Afc. In this model fi < 0.267 and Afc w 0.2. Therefore the condition 
is not satisfied. 

However, the plots in Ref.jBJ show steep increase of power for small k. The expected flat spectrum, 
which does not mean scale invariance, is not seen. This is no contradiction. When they plot the figures, 
they plot only discrete numbers of k modes and the plots do not have resolution less than Afc. This 
can be seen by looking at the shape of the resonance peak, which does not show the expected round 
shape. This means that the resolution is less than Afc. Since the mode k « Afc is inside the Hubble 
scale after the resonance becomes efficient, actually only one mode outside the horizon scale is plotted in 
those figures. Therefore it is not strange that the expected plateau is absent in their plots. As mentioned 
before, this n — 4 spectrum is common in the situation that there are random fluctuations of a large 
amplitude on small scales and seems to refute the claims of 0. 
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